PROCEEDINGS 

°F SCIENCE 



o 

O 
(N 

> 

o 
O 

(N 



00 
(N 



O 



Scaling and low energy constants in lattice QCD 
with A^f = 2 maximally twisted Wilson quarks 



Petros Dimopoulos, Roberto Frezzotti* and Gregorio Herdoiza* 

Univ. and INFN of Rome Tor Vergata, Via delta Ricerca Scientifica 1, 1-00133, Rome, Italy 
E-mail: 



{ dimopoulos , frezzotti, herdoiza} @roma2 . inf n . it 



Carsten Urbach 

Theoretical Physics Division, Dept. of Mathematical Sciences, University of Liverpool 
Liverpool L69 7ZL, UK 

t-mai I: Icarst en ■ urbach@physik . hu-berlin . de| 



Urs Wenger 

Institute for Theoretical Physics, ETH Ziirich, CH-8093 Ziirich, Switzerland 



E-mail: wenger @itp . phys . ethz . ch 



for the ETM Collaboration 

We report on the scaling of basic hadronic observables in lattice QCD with Nf ~2 maximally 
twisted Wilson dynamical quarks. We give preliminary results for some of the Gasser-Leutwyler 
low energy constants, the chiral condensate and the average mass of u and d quarks. 



The XXV International Symposium on Lattice Field Theory 
July 30 - August 4 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Scaling and low energy constants in N{ = 2 Mtm-LQCD 



Roberto Frezzotti and Gregorio Herdoiza 



1. Introduction and setup 

In this contribution we report on the scaling of basic hadronic observables and present prelim- 
inary results for some of the Gasser-Leutwyler low energy constants in lattice QCD with Nf = 2 
dynamical quarks. We choose the lattice formulation with tree-level improved gauge action and 
maximally twisted Wilson quarks which can be efficiently studied by means of state-of-the-art 
simulation algorithms, such as the one we adopted and leads to physical observables free of 
0(a) cutoff effects In this way it is possible to study the theory with pion masses down to about 
300 MeV, three different lattice resolutions and spatial lattice sizes of 2-3 fm. The main param- 
eters of the simulations employed for the present analysis are summarized in table 1 of Ref. [Q]. 
For all ensembles we have wipsL > 3, with the lowest values (namely 3.0, 3.3, 3.3 and 3.5) being 
obtained in the ensembles C6, Bl, CI and C5, respectively. For details about our lattice setup, the 
evaluation of quark propagators and any undefined notations we refer to Refs. ||l|, 0]. 

1.1 l\ining to maximal twist 

The values of K in table 1 of Ref. [Q] result from implementing maximal twist as discussed 
here. The formal definition of maximal twist for Wilson quarks reads: wtr = and /Ir = 0(a*') for 
all lattice spacings a as a —^0, with renormalized mass parameters 

Hr = Z^/I = Zp V 1 "JR = Z^o {m - merit) = ZAZp^mY>CAC ■ (1-1) 

At the non-perturbative level any legitimate estimate of the critical mass, niait, is affected by terms 
of 0{ah^Q^) and 0(a/iAQCD)> which however do not invalidate the definition of maximal twist. 
Following Refs. ||5|, ^, for each j8 (and /i) one can define maximal twist by demanding ^ 



a'MpcAc(j8,At) = 



„ =0, (1.2) 



at values of xq and L so large that the (charged) one-pion state dominates the correlators on the 
r.h.s. This prescription fixes the lattice artifact 0{aA^Q^) in the critical mass in such a way that, 
provided /i > a^Agf^p, the dominating (relative) cutoff effects left-over in physical observables ^ 
are expected to be numerically as small as 0{a^A^Q^). A detailed analysis ^ shows that such 
cutoff effects are actually products of (two) terms of order ^Aqcd. «M (negligible for small /i) or 
a^Agf-.p/i^^ (higher order in a, but enhanced for small ji). 

In practice, as we are interested in simulations with pL > /Ilow> in order to minimise the work 
for the tuning of K, we choose to impose the condition ( pr2| ) only for }JL = /Ilow- With this choice, in 
the region where /i < Aqcd we expect [^, |I|] the numerically dominating cutoff effects on physical 
observables to be modulated by factors of /^low //^ as /X is varied (for instance one can have contri- 
butions from terms of order (aAqcdMlow/m)^)- With K fixed according to the criterion described 
above, maximal twist is implemented properly, provided 

M ^ Row > Cfl^ Aqcd > (1-3) 



'Here quark bilinears are written in the unphysical quark basis (x, X) where the Wilson term has its standard form. 

^Concerning the ;r''-mass, in Ref. |Q| it is argued that the possibly large O(fl^) artifact on this observable is merely 
due to the large value taken by a (continuum) matrix element present in the Symanzik expansion of all lattice correlators 
where the ;r" state contributes, rather than to the presence of dimension six operators with large coefficients in the 
Symanzik effective action. This implies that this cutoff effect represents an exceptional, though important case. 
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where the coefficient C has to be learned from numerical experiment. In our Nf = 2 setup, C ~ 2 
(assuming Aqcd = 250 MeV) is essentially determined by the condition that MC simulations ex- 
hibit no metastabilities (see e.g. Ref. [|8|], sect. 5.5). At fixed /Ir a smooth approach to the continuum 
limit is to be expected if maximal twist has been realized at (approximatively) the same value of 
Mr low in physical units for all the considered lattice resolutions. 

Our practical implementation of maximal twist is illustrated by Fig. |l|a, where we show 
mRro mpcAC vs. pL^rQ oc /i for j8 = 4.05, 3.9 and 3.8 (use of renormalized quantities eases the 
comparison). Details on ro/a and Zp (actually Zp(MS,2 GeV)) are given below. 

For j8 = 4.05 and 3.9 we could fulfill our criterion for maximal twist at /Irlow — 0.047 
with good statistical precision: as shown in Fig. [I]a, mR//iRLOw = -ZkmpcAc/MLOW is consistent 
with zero within the statistical error, that we call e / /Ilow- As a rule of thumb, we demand e / jJ-how 
to be such that, numerically, aAgcD^/MLOW ^ 0.01. Considering the form of the d = 5 term 
in the Symanzik effective Lagrangian one finds in fact that aAqcd^/M is the expected order of 
magnitude of the unwanted (relative) cutoff effects that may contaminate physical observables 
if, due to numerical error, mpcAc(j8,A^) takes the value eZ^\ rather than zero. Having no a priori 
control on the coefficients of this order of magnitude estimate, we can only learn from the numerical 
experience of a scaling test (involving our statistically most precise observables, i.e. fys and mps) 
whether our rule of thumb yields a sufficiently accurate tuning to maximal twist. 
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Figure 1: (a) mRro = Z^Z^ 'mpcAc?-o vs. fi^ro and (b) ro(;U)/ro vs. (^ro)^, with ro = lim^^o '"o(m) (see 
text). In both plots data for )3 = 4.05, 3.9, 3.8 and L ~ 2.1,2.1,2.4 fm (respectively) are shown. 



At j8 = 3.8 our implementation of maximal twist was not as precise as wished due to the 
long-range statistical fluctuations that, for < O.Itq \ were observed in our MC simulations (see 
Ref. [0] for information on autocorrelation times). These fluctuations, which increase in MC-time 
length and (weakly) in amplitude as /i decreases, have a major impact on those observables, such 
as the plaquette and ampcAC, that are not continuous at the Singleton-Sharpe (I**' order) lattice 
phase transition [||, ^], rendering somewhat problematic the estimate of their statistical errors. This 
is manifest in Fig. |l]a from the two non-coinciding points obtained for the lowest /Ir at j8 = 3.8, 
coming from two independent simulations at slightly different L (ensembles A I and A5 in tab. I 
of Ref. iQ]). In this case, the quoted (possibly underestimated) statistical errors are larger than 
required by our rule of thumb. 
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1.2 Evaluation of ro 



In our scaling analyses we employ the Sommer scale ro 1 10] to eliminate the lattice spacing a. 
This is meant only as an intermediate step. In the end ro will be eliminated in favour of {i.e. fys 
at the physical point) to get the scale for all dimensionful quantities. The value of ro for the various 
ensembles was obtained with a better than 0.5% accuracy, starting from Wilson loops made out of 



HYP-smeared temporal links [ ]11[ ] and APE-smeared spatial links. Employing several interpolating 
operators for static quark-antiquaik (QQ) states, corresponding to different spatial smearings, leads 
to a matrix of correlators, from which the (lowest) levels of the potential VQg(r) are estimated 
by solving a generahzed eigenvalue problem. A fit to the r-dependence of the ground state gives 
ro(/x)/a. For each j8, the latter shows a rather mild -dependence, which is well described by a first 
order polynomial in /i^. The following values, extrapolated to the chiral limit and denoted simply 
as ro/a, will be used in this work (those at j3 = 4.05 and j8 = 3.8 must be viewed as preliminary): 

''o/a|/3=4.05 =6.61(3), ro/a|^=3.9 = 5.22(2), ro/a||3=3.8 = 4.46(3) . 
Note that a fit by a second order polynomial in jj. gives compatible results. An overview of the 
results for ro(/x) for all jS's is given in Fig. |l|b, where the axes are normalized in terms of the 
appropriate chirally extrapolated ro. For j3 = 3.8 the point at the lowest /i-value (ensemble Al) 
is preliminary and carries a still poorly estimated statistical error. It was hence not used in the 
extrapolation to the chiral limit. More details will be given in a forthcoming publication. 



2. Charged pion sector: general remarks 

In the charged pseudoscalar (PS) meson sector our raw results for fps and nip^ from simula- 
tions at j3 = 4.05 and 3.9 with L > 2 fm, once expressed in units of ro, exhibit an excellent scaling 
behaviour, see Fig. 4 of Ref. [Q]. Also the corresponding results at j3 = 3.8, in spite of the un- 
certainties on the implementation of maximal twist and the estimate of statistical errors, appear 
consistent with a very good scaling behaviour in the mass range mps = (350 600) MeV. These 
findings are in agreement with the expectation that cutoff effects are small in the absence of O(fl') 
artifacts. In particular in the charged PS-meson sector it is known ||5|, ^ that, for small quark 
masses /Xr, mp^ differs from its continuum counterpart only by terms of 0(a^/iR) and O(a^), while 
/ps has discretization errors of 0{a^). This property holds for any volume (sufficiently large to 
make pions much lighter than other states) and is not affected by the lattice artifact on the neutral 
PS-meson mass. ^ In other words, to order the cutoff effects on rtip^ and /ps are like in a chirally 
invariant lattice formulation. The use of continuum chiral formulae to describe the volume and 
quark mass dependences of our data in the charged PS-meson sector is thus well justified. 

2.1 Estimates of Zp and of the renormalized quark mass 

In the study of the scaling of the renormalized quark mass /Ir one needs the renormalization 
constant Zp{^;aq) at a common scale q for all the considered values of j8. This renormalization 
constant, as well as the scale-independent one Z^(j8), was also employed in Fig. [l]a. ^ Preliminary 

The relation of a possibly large and negative lattice artifact on Q to metastabilities in MC simulations has been 
discussed in Refs. [jlj ^]. Such metastabilities are not observed in the simulations we consider here [^. 

^For the renormalization constants we keep the names they are given in the literature for the standard (untwisted) 
Wilson quark lattice formulation, as obviously their values do not change with respect to the untwisted Wilson case. 
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0(a) improved estimates of the renormalization constants of quark bilinear operators are reported 
in Ref. [|l2|]. There the results for Za are rather precise (at the 1.5% level), while the quoted 
uncertainties on Zp are still substantially larger. 

In view of this situation, in our scaling analysis of }JLr we do not use the values of Zp quoted 



in Ref. [12]. We rather employ the scale- and scheme-independent ratios Zp(j3;<3^)/Zp(j8ref;fl'ref^), 
which we extract with a statistical accuracy of ~ 1% from the relation (exact up to 0{a^) terms) 

Zp{l3;aq) /Zp{l3rcf;ars{q) = Ai(j8;mpsro = l;L/ro ~ 5) //x(j8ref;mpsro = l;L/ro ~ 5) , (2.1) 

in order to compare the values of /Xs(j8;a^)Zp(j3ref;aref?) at different values of j8 and of the PS- 



meson mass. In eq. (2.1), a/x(j8;mpsro = l;L/ro ~ 5) is the value of ajx for which, at a given j8 
and L/ro ~ 5, one finds wipsro = 1. ^ The scaling is obviously not affected by the extra overall 
factor Zp (j8ref;aref^). which is removed in the end (in order to give an idea of the values of the 
renormalized quark masses). In the following we choose jSref = 3.9, where we most reliably know 
(see Ref. [|l^]) Zp(j8ref ; a^eiq), we evaluate for j8 = 4.05 and j8 = 3.8 the Zp-ratios from eq. ( ^!11 ) and 
we finally obtain estimates of pLR{^;aq) atq = 2 GeV in the MS scheme for all jS's. An important 
drawback of this method is the fact that cutoff effects stemming from data obtained at different 
j8 -values mix up in the quark mass renormalization, which in general may fake the genuine a- 
dependence and cast doubts on the reliability of any continuum extrapolation. Nevertheless, since 
our data for j3 = 4.05 and 3.9 show no statistically significant cutoff effects in the relation of /ps to 
and in the values of /Xp(j8;a^)Zp(j8ref;aref^), we can obtain an estimate of the continuum hmit 
of Hp (in units of ro) by taking an average of the results at these two j8 -values. A similar remark 
holds for the determination of Bq and the chiral condensate in sect. Q As discussed below, for all 
observables we associate to our continuum limit estimates a conservative systematic eiTor obtained 
by comparing them to the corresponding results from data at jS = 3.8. 



2.2 About taking the continuum, tliermodynamical and chiral limits 

It is well known that matching simulation data obtained at finite L and for mps > 300 MeV 
to the physical pion point requires a delicate analysis. The impact of residual O(a^) effects (even 
if small) on such an analysis might be enhanced if the continuum limit is performed as the last 
step. It is therefore advisable to perform first a continuum limit extrapolation at different (suitably 
chosen) fixed physical conditions and then use Chiral Perturbation Theory (;tPT) to correct for 



finite size effects and reach the physical pion point. This is the strategy we follow in sect. 4.1.1 



based on the results of the scaling test presented in sect. 3.1. As far as one is concerned with data, 
such as those at j3 = 4.05 and 3.9, where no significant cutoff effects are observed, a conceptually 
equivalent approach is that of performing a combined analysis of all data by means of continuum 
XVY formulae. The outcome of this approach, which is detailed in Ref. [Q], is summarised in 
sect. For comparison we also discuss in sect. [4. 1.3| the results of fits to continuum ;^PT 

formulae where the data corresponding to different lattice spacings are treated separately. 



^This value of mpsro was chosen so as to lie in the region where errors from interpolation to the reference mass and 
cutoff effects are smallest. 
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3. Charged pion sector: scaling test 

Here we analyse the scaling behaviour of the charged PS -meson decay constant and the renor- 
malized quark mass as a ^ at fixed values of mpsro and L/r^. The renormalization and scaling 
conditions are as follows: (i) rQFgg(ro) = 1.65 (with Fqq the static interquark force), which allows 
to trade = 6/j8 with ro/a ; (ii) mpsro = constant G {0.7,0.8,0.9, 1.0, 1.1, 1.25}, which eliminates 
/i in favour of mps ; (iii) fixed spatial volume L^^^ ~ (2.2 fm)^, which corresponds to L/ro ~ 5. 

3.1 Scaling and preliminary continuum limit estimates 

The typical relative statistical errors are (conservatively) estimated to be about 0.7% for mpsro, 
1.0% for /psJ'o and 1.5% for pL^VQ. In the latter case, besides the error on ro/a, only the uncertainty 



(typically about 1.2%) on the Zp-ratios (see eq. (2.1)), is taken into account and hence shown in 



Fig. ^. The error on Zp{^sf',civefq) is omitted here, as it plays no role, while in sect. ^ it is taken 
into account, as it matters for the final values of Bq, the chiral condensate and niud- 

The condition (i) is immediately fulfilled once all quantities are expressed in units of tq. As 
for the other two conditions, we first implement the condition (iii) by "moving" via resummed %¥! 
formulae (see Ref. [14|) the data for /ps and mps from the simulation volume to the reference 
one L^gj. As simulations were done with L G (2.0,2.4) fm, the numerical change in the data is 
very tiny (never larger than 0.7%) and thus statistically almost irrelevant. Then, in order to match 
the reference values of mpsro given in (ii), we perform interpolations (in few cases also short 
extrapolations) of the values of /ps^o and /iro. For this purpose we try both (low order) polynomial 
and ;^PT-inspired fits to the data. The spread among the different fits with good whenever 
statistically significant, is added linearly to the interpolation error. In this way one ends up with the 
blue filled circle data points (and errorbars) in the Figs. § and |3[ 

For each value of mpsro we obtain preliminary estimates of the continuum limit values of /ps^'o 
and /iRro by fitting to a constant (red line in the figures) only the data points from simulations at 
j8 = 4.05 and j3 = 3.9. The results, with only statistical error from the continuum extrapolation, 
are shown in red in Figs. |2| and ^. In all cases the difference between the result of the continuum 
extrapolation and the central value (blue filled circle point) for j8 = 3.8 is taken as an estimate of 
the systematic error (indicated with a green cross placed at a slightly negative value of a) on the 
continuum limit result. Since the latter is obtained by simply taking a weighted average of the 
results at j8 = 4.05 and 3.9, the introduction of such an error appears necessary. Moreover, the 
quality of the data in the figures suggests that the way we estimate this systematic error is rather 
conservative. Finally, the green crosses (and errorbars) appearing in Figs. ^ and ^ at (a/ro)^ ~ 0.05 
show, whenever the displacement is larger than one standard deviation, where the actual data would 



move if one were correcting for the leading effect of the numerical error (denoted by e in sect. |1.1| ) 



in imposing nipcAC = 0. We refer to the results of such a correction (explained in sect. |3.2D as to 
"data moved to maximal twist". 



3.2 About "moving data to maximal twist" 

The procedure of "moving data to maximal twist" can be seen as a way of testing whether 
statistical errors in the tuning to maximal twist have any significant impact on the observables of 
interest. If one has Z^mpcAC = £ (rather than zero), the effective twist angle in the Symanzik 



6 



Scaling and low energy constants in N{ = 2 Mtm-LQCD 



Roberto Frezzotti and Gregorio Herdoiza 



mp5rQ = (0.8. 1.0, 1.25) 
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,2 



0.04 0.05 0.06 



(a) (b) 
Figure 2: Scaling plots of /psro vs. {a/ro)^ at fixed values of mpsro (increasing from bottom to top). See 
text for a detailed explanation of the symbols. 



Lagrangian is a = 7i/2 — 6 with land = = e/jj.. Treating e as an 0(a°) quantity and 

neglecting 0{ad) and O(fl^) effects, the consequences of a deviation from maximal twist are in 
general twofold. First, the effective renormalized quark mass becomes Mr = (/i^ + m|)^/^ = 
Ha/cosd. Second, in all operator matrix elements any operator formally non-invariant under 
axial-T^ transformations must be reinterpreted consistently with the twist angle being a = n/2— 6 
(rather than 7i/2). As a consequence, once the hadronic states of interest are correctly identified, 
simple 6-dependent formulae can be derived that allow to extract properly the matrix elements 
of operators non-invariant under axial-T^ transformations. For instance, the (renormalized) op- 
erator Zyxy^'^^X (written in the basis where the Wilson term has its standard form) coincides 
with cos ^ + sin dV^jf, where ^ (^;?r) is the physical axial (vector) current. Consequently 
one finds < n^\ZvX7^i"^^X\^ > = cos0 < 7r'|A^^|n > |m« +0{ad,a^) , which implies 

mps/psk^ =< 7l^\ZvZ7o^^Z\^ > \m^.nJcose + 0{ad,a^). 

From these arguments it follows that nips-data need not be moved, while the data for ^LiR and 
fps are "moved to maximal twist", up to 0{ad,a^), by dividing them by cos d, where 6 is obtained 
from the (interpolated) actual values of e and n in our data sets. It turns out that the change implied 
by this correction is larger than one standard deviation (of the uncorrected data) only for the two 
smallest reference values of mps^o at j8 =3.8. After such correction all the data points at j8 = 3.8 
(including the two most chiral ones) fit very nicely into the scaling picture suggested from the 
corresponding data at j8 = 3.9 and j8 = 4.05. 

However, any data where the above correction was statistically relevant can be hardly used 
for a continuum extrapolation, because the residual (uncorrected) 0{ad) lattice artifacts depend on 
the details of the numerical errors at the different j8 -values and hence need not to scale as a ^ 0. 
Nevertheless, at a more qualitative level, the procedure of "moving data to maximal twist" confirms 
that, even at j3 = 3.8 (i.e. a ~ 0. 1 fm), the lattice artifacts in the charged PS-meson sector appear 
to be quite small, provided maximal twist is implemented precisely, e.g. according to the criteria 
of sect. |1.1| . In fact such a precise implementation is rather hard to achieve in the case of our 



7 



Scaling and low energy constants in N{ = 2 Mtm-LQCD 



Roberto Frezzotti and Gregorio Herdoiza 



m r = (0.7, 0.8,0.9) 



m„,r = ( 1.1, 1.25) 



0.01 0.02 0.0.^ 0.()4 0.05 0.( 




(a) (b) 
Figure 3: Scaling plot of /^Rro vs. (a/ro)^ at fixed values of mpgro (increasing from bottom to top). See text 
for a detailed explanation of the symbols. 



simulations at j8 = 3.8, as we found Tint(fl'"^PCAc) = O(IOO), see sect. 1.1 and Ref. 



4. Charged pion sector: physical results 

Here we discuss the description of our mps- and /ps-data by means of continuum xVY for 
Nf = 2 QCD and give preliminary estimates of the low energy constants (LEC) Is, U, Bo and /o, ^ 
as well as of the average light quark mass, m„j, and the chiral condensate. 



We employ the following continuum NLO ;^PT formulae [13, 14] to simultaneously describe 



the dependence of mps and fps on the bare quark mass (/i) and on the finite spatial size (L): 



m 



2 

PS 



(L) = 2Bo^Kl{L) [1 + ^ ln(2BoAt/Al)] , (4-1) 



/ps (L) = foKf{L) [1 - IE, ln(2BoAt/A4)] , (4-2) 

where t, = IBojX / {AnfoY ^^^^ Kmj{L) account for finite size (FS) effects. The LEC are related 
to the parameters A3 4 introduced in eqs. (|4. l[)-( |4!2| ) through /a 4 = log (A3 4/ m^± ) . The expressions 



foxKmj{L) at NLO (denoted asGL)read /<:,^^(L) = \ + {t,gi{X) wdKy^iL) = l-2t,gi{X), 



where A = \/2Bq}JlL^ and gi{X) is a known function. A convenient way to include higher order 



XVI terms in the description of FS effects on mps and /ps is provided by the formulae of Ref. [14] 
(denoted by CDH). In the following we always use the CDH expressions for K,„ f{L), which turn 
out [Q] to describe well our data at different volumes (and for mps better than GL formulae). 

In all the analyses below physical units are introduced as follows: the experimental values, 
fji = 130.7 MeV and ntj^ = 135.0 MeV, ^ are exploited to determine first the /x-value, cor- 
responding to the "physical point" through nips//ps|^;i = nij^o/fji, and then (depending on the 
analysis) the value of rofj^, or afj^. The latter allows to express all quantities in units of fj^. 



^We use the convention /q — ViFq, i.e. the normalization fjc = 130.7 MeV. 

^The ;r"-mass input is chosen in view of the absence of electromagnetic effects in our lattice QCD simulations. 
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4.1 Chiral fits and continuum estimates of LEC 

As anticipated in sect. [2.2[, given the good scaling behaviour of our data, we present prelim- 



inary continuum estimates of LEC and rriud that stem from two analyses - see sects. 4.1.1 and 
|4.1.2| . In sect. [4.1. 3| we give the corresponding results that are obtained at two fixed lattice spac- 
ings, = 4.05) and a(j3 = 3.9). The latter results will serve mainly for estimating (in the way 
detailed below) the systematic error due to residual O(a^) cutoff effects. 

Concerning other systematic errors, we note: (i) residual uncertainties in the CDH-formulae 
for FS effects are small, compared to other systematic errors - see the discussion in Ref . [Q] ; (ii) the 



possible impact of NNLO corrections to the formulae (|4. 1D-([4.2|) for the quark mass dependence is 
minimized by taking out of the analyses the points with highest values of mps (we checked that fits 
are stable if we leave out data with mps > 500 MeV - more details in a forthcoming publication). 

4.1.1 ;^PT fits in the continuum 

Following the strategy presented in sect. one can estimate the continuum limit values of 



/psro and jXurQ at fixed values of m^^rQ and L^^f/rQ, as reported in sect. Considering /ps?"o and 
mps^o as functions of Huro allows for direct use of the ;^PT formulae (|4~Tl)-(p^ to bring our data 
from Lref = 2.2 fm to infinite volume and parameterise their quark mass dependence. Leaving out 
the point corresponding to mpsn) = 1 -25, the data are well described by our fit ansatz and we obtain 
the following values for the fit parameters: 

IBoro = 12.0(3)(7), h = 3.67(12)(35) , 
Mo = 0.266(3)(10), k = 4.69(4)(11), (4.3) 

where the renormalized quantity Bq = ZpBo is given in the MS-scheme at the scale q = 2 GeV (see 
sect. 2.1). The ;^^/dof of the fit is 0.28. Using the experimental input at the physical pion point we 



find ro = 0.433(5)(16) fm. Inserting the values of ro/a (see sect. |1.2[ ), one gets the estimates 

a\p=4.05 = 0.0655(8)(24) fm, a\p^^ g = 0.0830(10)(31) fm, a|/3=3 g = 0.0970(13) (37) fm. 

In the above results, except for the case of 2Bo'"o> the second error comes entirely from the 
systematic uncertainty due to residual O(a^) cutoff effects: for each quantity this error is conser- 
vatively taken as the maximum of (i) the uncertainty resulting from the propagation through the 



chiral analysis of the systematic error associated to the continuum estimates as derived in sect. |3J 



and (ii) the spread of the results obtained separately at j8 = 3.9 and j8 = 4.05 (by fitting to the same 



;fPT ansatz as here - see sect. 4.1.3). In practice the maximum is usually given by the former of 



these two systematic error estimates. Note however that for Bq, niud and the chiral condensate, an 



additional systematic uncertainty coming from Zp (as quoted in Ref. [|12p) is added in quadrature. 

As a check, we also study the decay constant /ps as a function of mps, i.e. with no reference 
to /Ir. The appropriate NLO ;^PT fit ansatz is obtained (ignoring NNLO corrections) by replac- 
ing 2B()/x with in eq. ( |4.2[ ). The resulting best fit parameters are /oro = 0.268(3)(12) and 



U = 4.82(4)(14). This is consistent with the values in eq. (4.3) and yields for the Sommer scale 



the estimate ro = 0.435(4)(15) fm. The second error again comes from (the propagation of) the 



systematic uncertainty in the "continuum extrapolation" of sect. 3. 1 



9 



Scaling and low energy constants in N{ = 2 Mtm-LQCD 



Roberto Frezzotti and Gregorio Herdoiza 





single-j3 fit 


single-/3 fit 


combined fit - see sect. [4.1.2 




3.9 


4.05 


3.9 4.05 


2aBQ 


4.85(4) 


3.87(6) 


4.87(4) 3.76(3) 


afo 


0.0526(4) 


0.0404(7) 


0.0527(4) 0.0411(4) 


A3//0 


6.36(26) 


7.20(48) 


6.41(26) 


A4//0 


11.59(19) 


11.81(31) 


11.51(21) 


zVdof 


7.8/6 


2.2/4 


12.0/12 


ajin 


0.00072(2) 


0.00054(2) 


0.00072(1) 0.00057(1) 


a [fm] 


0.0854(6) 


0.0656(10) 


0.0855(5)(31) 0.0667(5)(24) 



Table 1: Results from fits at single j3 -values and from the global combined fit of sect. 4.1.2. 



4.1.2 ;^PT analysis combining j8 = 3.9 and 4.05 



As argued in sect. given the absence of statistically relevant cutoff effects, one can com- 
bine the data (in lattice units) for nips and /ps coming from our simulations at j8 = 3.9 and j8 = 4.05 
and perform a global combined fit based again on eqs. (|4. 1[)-( |4!2| ). Such a combined fit has six free 
parameters which can be taken as (aBo)|/3=3.9> (afio)|/3=4.05> («/o)|/5=3.9> (a/o)|/3=4.05> A.3//0 
and A4//0. Details and plots concerning this analysis can be found in Ref. [Q]. The outcome is 
summarised in table |l|. As for the estimates of a, the second error we quote is an estimate (see 
sect. [4.1.1| ) of the systematic uncertainty due to residual cutoff effects. We note that, leaving out 
the cases for which mps > 500 MeV, nine ensembles of gauge configurations enter this analysis: 
five at j3 = 3.9 (Bi to B4 and Be) and four at jS = 4.05 (Ci to C3 and C5). 

4.1.3 Independent ;^PT analyses at j3 = 3.9 and 4.05 and comparison 



These analyses, always based on eqs. (4.1)-(4.2), are closely analogous to that presented in 



Ref. [jl|]. The free parameters can be taken to be qBq, afo, A3//0 and A4//0, the best fit values 
of which are given in table [I] and illustrated in Fig. Q For the case of j8 = 3.9 the only difference 
with respect to Ref. [|l]] is the increased statistics and in particular the presence of the ensemble 
Be- For both j8 = 4.05 and j8 = 3.9 the same ensembles are considered as in the case of the global 
combined fit of sect. [4.1.2| . It should be noted that in the analysis of sect. [4.1. 1| the ensembles C5 and 
Be were not employed, in order to avoid to heavily rely on the CDH-formulae in the "continuum 



extrapolations" discussed in sect. [3.1[ On the other hand using these two ensembles, in particular 
Be (which corresponds to L ~ 2.7 fm) increases the statistical information. 

In table ^results from the independent analyses at j8 = 3.9 and 4.05, the combined fit of 
sect. [1.1. 2[ and the fit in the continuum of sect. [1.1. 1| are compared (the second error, whenever 
quoted, is the systematic one previously discussed). The agreement between the last two columns 
of table ^ is good, even within statistical errors only. We recall that the ensembles Be and C5 do not 
contribute to the results in the last column. Were these ensembles not used also in the combined fit 
of sect. [4.1. 2[ , the results in the third column of table § would get, as expected, even closer to those 

'*From this analysis, without using ro, we find a\p^2 (j/a\p^4 Q^ = 1.28(1), which is quite close to 
(ro/fl)|^=4 05/(''o/fl)|j3=3.9 = 1.27(1). Moreover, the comparison of the ratios /o/6o|j3=3.9 and /o/So|j3=4.o5 provides 
a rather precise estimate of the ratio Zp{fi — 3.9;ap^^gq)/Zp{fi — 4.05; 0^3=4 os?)- 
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(a) (b) 
Figure 4: Independent jPT fits at j3 = 4.05 and ji = 3.9: FS-corrected data a la CDH for (a) rQUip^/pi^ 
and (b) ro/ps vs. /iRro are shown, with the corresponding best fit curves to eqs. (4.1)-(4.2). The points at 
^RrQ^O.ll were not included in the fit. 





single-fit at j3 = 3.9 


single-fit at j3 = 4.05 


combined fit 


fit in the continuum 


h 


3.41(9) 


3.66(14) 


3.44(8)(35) 


3.67(12)(35) 


u 


4.62(4) 


4.66(7) 


4.61(4)(11) 


4.69(4)(11) 





11.6(3)(5) 


11.9(3)(5)-12.1(3)(5) 


11.6(3)(7) 


12.0(3)(7) 


fon) 


0.275(2) 


0.267(5) 


0.273(3)(10) 


0.266(3)(10) 


ro [fm] 


0.446(4) 


0.434(7) 


0.444(4)(16) 


0.433(5)(16) 



Table 2: Estimates of LEC and ro from the analyses of sects. 4.1.3 (1*^' and 2"*^ column), 4.1.2 (3"* column) 
and 4.1.1 (last column). The values of 2roBo are obtained using Zp at j3 = 3.9 JT^ ] and Z^-ratios from the 
analyses of sects. 4.1.2 and 4.1.1: at j3 = 4.05 we quote the two (similar) results obtained in this way. 



in the last column: h = 3.66(9), U = 4.64(5), 2roB() = 1 1.9(3), /oro = 0.270(3), ro = 0.439(5) fm. 



For reviews on determinations of the LEC we refer to Refs. [113, 16]. 



4.2 Continuum estimates of ruud and the chiral condensate 

Our preliminary results in the MS-scheme at scale q = 2 GeV are obtained using the RI-MOM 



estimate [ ]12[ ] of Zp at j8 = 3.9 (converted to MS by using NNNLO perturbation theory) and the 



Zp-ratios from the analysis of sect. [4.1.2 (combined fit) or [4.1.1| (fit in the continuum). The average 



light quark mass mj,^(MS,2GeV) from the fit in the continuum reads = 3.43(9) (23) MeV. 
From the combined fit we get niud = 3.62(10) (23) MeV (using all the ensembles) and niud = 
3.49(13)(23) MeV (when excluding ensembles Be and Cj, i.e. with the same set of ensembles 
as for the fit in the continuum). When considering data at j3 = 3.9 only (as in sect. [4.1.3| ), we 
obtain niud = 3.62(13) (23) MeV This value is compatible with the one coming from the partially 
quenched analysis of Ref. niud = 3. 85 (12) (40) MeV. 

The chiral quark condensate is obtained through (^g)(MS,2GeV) = -(l/2)/(^Bo(MS,2GeV). 
From the fit after the continuum extrapolation we get \ {qq)\^^'^ = 272(4) (7) MeV, while from the 
combined fit analysis we obtain \ {qq)\^/^ = 267(4)(7)MeV (and \ {qq)\'^l^ = 270(4)(7)MeV when 
the ensembles and C5 are not used). These results are in agreement with an independent estimate 



at j3 = 3.9 from the e regime - see Ref. [18]. 
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Figure 5: myro and nij^ro vs. (mpsro)^ for different lattice resolutions and L in the range 2. 1 2.4 fm. 



5. Scaling analysis of other hadronic observables 

Here we briefly report on the scaling behaviour of other hadronic observables. Even if at 
jS = 3.8 preliminary data are available only for the meson vector mass and the physical volumes at 
the various lattice resolutions (L ~ 2.1 fm for j8 = 4.05,3.9 and L ~ 2.4 fm for j8 = 3.8) are only 
approximatively matched, this overview suggests that when employing maximally twisted Wilson 
quarks the cutoff effects on physical quantities are in general as small as expected on the basis of 
0(a) improvement. 



5.1 Nucleon and vector meson masses 

In Fig. ^, we show our data (in units of ro) for the vector meson "mass", my, and the nucleon 
mass, /TIN. as a function of nip^. For the former one (see Ref. [ [I9| ] for more details), the scaling 
appears to be good, although within statistical errors that become of few percents when nips ~ 
300 350 MeV. In this case, we remark that what we call the "mass" may differ from the actual 
vector meson mass (the one that becomes nip as mps — > rrijc), owing to the effect of virtual p-nn 
mixing. ^ Due to the finite volume, the decay of the p -meson (at rest) into real Tin states is instead 
forbidden, even when 2mps < my (which roughly speaking happens in the region (mpsro)^ < 1). 



The data for the nucleon mass are also prehminary (see Ref. [20]): at j8 = 4.05 only two data points 
are available, which however appear to agree quite well with the data points at j8 = 3.9. 



5.2 Charmed observables in partially quenched setup 

To compute charmed observables we adopt a partially quenched (PQ) setup, which is detailed 
in Ref. [21]. With the notation of Ref. [ |22| ] for valence quarks, the correlators for D (cd) and 
Ds {cs) charged mesons are computed with Wilson parameters = = — rc = 1, in order to have 



^The correction is estimated [jl^ to be almost irrelevant within our present statistical errors: at our most chiral 
point, where mps ~ 300 MeV, it amounts to ray — "Jy ~ 0.05my. 

"'The proof of automatic 0(a) improvement given in Ref. | |2^ ] for a fully unquenched mixed action framework goes 
through also in the PQ mixed action setup of Ref. where the quark masses of the valence s and c quarks are taken 
finite and hence different from those (infinite in the present case) of the corresponding sea quarks. 
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Figure 6: foro, /ds'"o (a) and m^jo (b) vs. moro for two lattice spacings (j3 — 4.05 and 3.9) and L^^IA fm. 



the same nice parametric scaling properties as in the charged pion sector. For the present consider- 
ations about scaling we choose close-to-realistic renormalization conditions: mpsro = 0.7092 and 
M.v/Mc = 0.082 (instead of the reahstic conditions nij^ro ~ 0.30 and IxJlJ-c — 0.088). In Figs. ^, 
mD varies with [Xg at fixed values of /x^/jUf and nipsro {i.e. }Xud)- The expected size of the domi- 
nating cutoff effects is {ajXcf ~ 0.1, since typically, a/if G [0.020, 0.033] and a^^ G [0.25, 0.40]. 
By comparing the preliminary results from j8 = 4.05 and jS = 3.9, we see that the observed scaling 
violations are not large: from 1-2% (hardly visible within statistical errors) for /d and mo, to 7-8% 
for /dj. Data from a third lattice spacing, j8 = 3.8, are currently being analysed and may allow for 
an estimate of continuum limit results. 

6. Conclusions and outlook 

We have reported on the scaling properties of lattice QCD with Nf = 2 maximally twisted 
Wilson quarks, in the unitary [Q] as well as in a partially quenched setup |[T^, Very good 
scaling properties are found in the light (charged) PS-meson sector and also for various charmed 
PS-meson observables, as well as in the vector meson and the parity-even nucleon channels, in 
agreement with the expectation of automatic 0(a) improvement of physical observables. Based on 
these findings, we presented ;^PT-based analyses suggesting that the determinations of the LEC's 
(among which those relevant for the chiral condensate) and of the average u,d quark mass that were 
obtained for one single lattice resolution in Refs. ^ |l^ are likely to be close to their continuum 
limit Nf = 2 QCD values. Moreover, combining the results presented here with those of Ref. [^], 
after appropriate chiral extrapolations, one finds an estimate of the ratio m^/ fji that is in good 
agreement with experiment. The framework of lattice QCD with maximally twisted Wilson quarks 
appears thus to offer good prospects for reliable computations of many physical QCD observables 
and weak matrix elements in the continuum limit. 
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